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O I Quantized nonlinear lattice models are considered for two different classes, boson and 

fermionic ones. The quantum discrete nonlinear Schrodinger model (DNLS) is our main 
objective, but its so called modified discrete nonlinear (MDNLS) version is also included, 
^ ■ together with the fermionic polaron (FP) model. 

^ , Based on the respective dynamical symmetries of the models, a method is put forward 

qh| which by use of the associated boson and spin coherent state vectors (CSV) and a factoriza- 

tion ansatz for the solution of the Schrodinger equation, leads to quasiclassical Hamiltonian 
equations of motion for the CSV parameters. The so obtained evolution equations are in- 
timately related to the respective evolution equations for the classical lattices, provided we 
account for the ordering rules (normal, symmetric) adopted for their quantization. 

Analysing the geometrical content of the factorization ansatz made for the state vectors 
invokes the study of the Riemannian and symplectic geometry of the CSV manifolds as 
generalized phase spaces. Next, we investigate analytically and numerically the behavior 
of mean values and uncertainties of some physically interesting observables as well as the 
modifications in the quantum regime of processes such as the discrete self trapping (DST), in 
terms of the Q-function and the distribution of excitation quanta of the lattice sites. Quan- 
tum DST in the symmetric ordering of lattice operators is found to be relatively enhanced 
with respect to the classical DST. Finally, the meaning of the factorization ansatz for the 
lattice wave function is explained in terms of disregarded quantum correlations, and as a 
quantitative figure of merit for that ansatz a correlation index is introduced. This index is 
given in terms of the norm of the difference between the true and factorized state vectors, 
and accounts for the quantum correlations of the lattice sites that develop during the time 
evolution of the systems. 
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1 Introduction 



Nonlinear lattice equations are abound in physical applications and their study constitutes 
a well developed field of research. In many cases where models are described by classical 
lattice differential equations this is done as a first approximation to what is essentially a 
quantum mechanical lattice system. Quantized nonlinear lattice equations are the main 
target of the most successful method of the quantum inverse scattering method (QISM) p5| . 
Recent attempts to treat quantum lattice equations with arbitrary number of lattice sites 
have led to such schemes as the so called number state method (NSM) |3^, |T^, and the 
Hartree approximation (see also for a discussion about QISM and NSM). On the 
other hand it seems highly desirable when treating quantum versions of classical lattices to 
be able to utilize in some sense the available classical solutions of the nonlinear evolution 
equations and to keep a conceptual framework wu.phase spaces, Poisson brackets, coordinate 
transformations etc. akin to that of the classical models. It is in that direction that this 
work focuses and provides an alternative to the above mentioned methods. 

Specifically we deal with lattices related to the nonintegrable generalized discrete non- 
linear Schrodinger equation (DNLS) (also known as the generalized discrete self-trapping 
(GDST) equation) [^] and with a modified version of it (MDNLS) ||21|. We also discuss a 
fermionic lattice associated to the fermionic polaron (FP) model [l^. The DNLS models 
quantized by the canonical quantization method lead to a configuration of / sites, described 
by tensoring the Weyl-Heisenberg (WH) algebra attached to each lattice site. The ordering 
rule adopted when the classical wave amplitude variables are substituted by boson opera- 
tors ^ is essential for the form of the obtained Hamiltonian operator. Our method of 
solving the resulting Schrodinger equation is based on the assumption that the lattice state 
vector can be factorized to a tensor product of coherent state vectors ([^, |17[ ^ ^ ^) each 
one living in the Hilbert space associated to each lattice site. This leads to a set of Hamilto- 
nian equations for the CSV coordinates, which provide an approximate quasiclassical solution 
for the exact quantum dynamics. As the CS are the most classical state vectors |3^, [l^ for 
boson systems the resulting equations are quasiclassical in nature and govern the motion 
of CS wavepackets. To elucidate the meaning of the ansatz we recall that the boson CS 
labels are coordinates of the canonical phase plane and therefore the quasiclassical dynamics 
occurs in the /-fold cartesian product of harmonic oscillator phase space. That dynamics is 
the classical dynamics if the normal ordering of the bosons is used and diverges from that 
for symmetric ordering. Therefore the CS mean values for normally ordered observables 
coincide with their classical values in the former case while their quantum mechanical fiuc- 
tuations are constant and minimum in the course of time evolution. To try the method to 
more general quantum systems we have also studied the FP model, which is first mapped by 
means of the Jordan- Wigner transformation to its equivalent XXZ model 0. The resulting 
spin Hamiltonian is expressed by tensoring the su{2) algebra generators attached to each 
site of the model. By the same token as in the boson case the quantum dynamics is studied 
by assuming factorization of the state vector of the system to a product of spin CS. How- 
ever the geometrical framework for the dynamics of the spin CSV's is different in this case. 
The su{2) CS space is a spherical surface with a noncanonical symplectic structure defined 
on it, as it results by looking at the so called ray metric of the CSV |^ |3^; the classical 
Hamiltonian equations generate a fiow in the /-fold cartesian product of spherical (or rather 
complex projective) coordinates. 
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The plan of the paper is as follows: in section two the general framework for the con- 
struction of CSV is outlined and the needed formulas for the special cases of the canonical 
and the spin CSV are provided. In the next section the DNLS model is quantized and the 
associated quasiclassical equations are derived and studied. Attention is paid to the prob- 
lem of discrete self trapping (DST) phenomenon in the quantum regime, by illustrating the 
situation with plots of the Poisson distribution and the Q-function of the quantized modes 
of the system. Next section explains the physical meaning of the factoring of the lattice 
state vector as adopted for all models, and gives a validity measure for that approximation. 
In the following two appendices we illustrate the same idea by taking up first the MDNLS 
quantum lattice and then the FP model. The final section contains a number of conclusions 
and offers some prospects. 

Throughout the text we consider h = 1 and we denote CS expectation values of any 
generator with brackets around it. 



2 Coherent state manifolds 

The notion of coherent state vector was essentially introduced in the early days of quantum 



mechanics by Schrodinger ||3^ , in his succesful attempt to construct a nonspreading Gaussian 
wavepacket for a quantum harmonic potential. The wavepacket center was evolving following 
the classical path while the dispersion was kept to a minimum compatible with the minimum 
uncertainty principle. The regeneration of the theory took place in the sixties and seventies 
when to a wealth of physical applications of CS a mathematical group theoretical foundation 
was also provided [T|, |2|, |^. 

For our needs here a brief introduction of the CS concept goes as follows: consider a Lie 
group Q, with a unitary irreducible representation T{g), g G Q, in a Hilbert space H. We 
select a reference vector |\I'o) G H, to be called the "vacuum" state vector, and let Qq C Q 
be its isotropy subgroup, i.e. for h E Qo^ T{h) |\l/o) = e^"^^^^ l^o)- The map from the factor 
group M. = Q / Qq to the Hilbert space ?i, introduced in the form of an orbit of the vacuum 
state under a factor group element, defines a CSV \x) = T{Q /Qq) |\&o) labelled by points 
X E M. oi the coherent state manifold. Coherent states form an (over) complete set of states, 
since by means of the Haar invariant measure of the group Q viz.dfi{x), x E A4, they 
provide a resolution of unity, 1 = djj,{x) \x){x\. As a consequence, any vector G Ti 
is analyzed in the CS basis, |\I') = jj^dii{x)'^{x) with coefficients '^{x) = (a;|\l/). 
What concerns us here mostly is the geometry of the CS manifold A4. Indeed by its very 
construction Ai inherits the structure of a Riemann manifold with in general non-constant 
curvature, which is also endowed with a complex structure of a Kahler manifold [^, namely 



it can be considered as a generalized phase space [0]. In the sequel we restrict ourselves 
to the case of a two dimensonal surface Ai for definiteness, although higher dimensional 
extensions of our statements are also possible Also we shall assume a generator G+ 

creating from the vacuum state the CS vector, i.e., 

|C)=Ar||C)=Arexp(CG+)|0) , (1) 

where C G C and J\f = (C||||C)~^ is the normalization factor. Occasionally we shall write 
G_ for the Hermitean conjugate of G+; they should both belong to the Lie algebra of the 
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group for which the CS is defined. Below we shall specialize to the cases of su{2) and the 
Weyl-Heisenberg groups but for the moment we proceed with the present general framework. 
The symplectic structure possessed by the state space is based on the existence of a 
canonical kinematical 1-form 6 = {(\d\() = ^{{G^)d( — {G^)d(). The derivation d = 
d^dC, + d^dC, acts on the state vectors, e.g., d |C) = (G+ — \ {G+)) |C) dC + (— § (G-)) |C) dC. 
The symplectic 2- form of is derived from the canonical 6 by derivation u = d6, and 
can be expressed in the form u = — (G-) {G+))d( A d(. 

Concerning the geometric features of our 2D phase space considered as a Riemannian surface 
with a distance function operating on it, we shall employ a meaningful metric tensor starting 
from the distance on the projective Hilbert space of rays \0, ( G C ||2^, 0. By \() we simply 
mean the set of CS e**^ \(), multiplied by an arbitrary phase factor. The choice of space of rays 
rather than the Hilbert space is in accordance with the quantum mechanical arbitrariness 
of phase of the state vectors. The finite distance V{\(i), 1(2)) between any two rays 
and IC2), which are associated with the normalized coherent states e**^^ and e*"^^ IC2) 
defined by 



IS 



^^^(ICi), IC2)) 



inf 



lCi)-e*^^|C2 



2|(Cl|C2)| 



(2) 



This is a proper distance function as it is positive definite and non-degenerate, it satisfies 
the triangular property and it is also gauge invariant. In addition its infinitesimal form gives 
the distance of two nearby coherent state vectors and provides the metric tensor on the 
CS-manifold HOI, i.e. 
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ids) 



ids) 



m\c+do, ic)) = i-i(c+rfcic)i 
(CI 1 1 10 5c (CI 1 1 10 



^c^dCIIIIC) 



d 



(CIIIIC) 

icr 



5icr 



(CIIIIC) 
^(CIIIIC) 



(CIIIIC) m 



(CIIIIC) 

dCdC = 9((dCdC 



dCdC 



(3) 
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The curvature scalar IT 



R = -%^5c5c(l^^cc) 



(5) 



following from the metric g involves higher correlations of the G± generators and is in general 
not constant. It is worth noticing that the basic geometric objects of both the symplectic 
and Riemannian structure endowed in Ai are given in terms of the so called symbols of the 
operators or their products 0, i.e. the coherent state mean value (C| G |C) = (G) of the 
corresponding operator G. Also the noncommutativity of the involved operators and their 
non-zero uncertainties in the coherent state basis is essential for the non-trivial geometric 
characteristics of the Ai manifold. In effect, the CS-manifold captures some genuine 
quantum mechanical features despite its classical character; this property of Ai is further 
manifested in the quantum mechanical evolution to be studied shortly. 
As the models to be studied in the next sections are of boson and fermionic/spin type with 
dynamical symmetries related to the WH and the su{2) algebras |T^, we now exemplify the 
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above construction for the CS of these algebras. For the boson algebra, 



[a, a^] = 1 [TV, a^] = [TV, a] = -a , (6) 

with = a'I'a the number operator, the vacuum state is the zero-photon state vector |0) — 
|0), and G+ = a^" the creation operator. This defines the boson CS 

oo n 

It is an (over) complete set of states with respect to the measure dix{a) = -e~l"l^(i^Q; for the 
non-normalized CS, and a & M. = WH/U{1) a; C is the CS manifold. Since a\a) — a \a), 
which implies the symbols (a) = a, (a)'^ = a, (N) = \af , Ai is the flat canonical phase 
plane with the standard line element ds"^ = dada. Also the symplectic 2-form uj = ida A da 
is associated to the canonical Poisson bracket {f,g} = i^dafdaQ — dafdag)- 
Next we come to the case of the su{2) algebra with commutation relations, 

[Jo,J±]=±J± [J+,J_] = 2Jo. (8) 

The vacuum state |0) = \j — j) is the extremal vector in the representation module of the 
algebra, with {\j m)}l^^_j, j = 1/2, 1, 3/2, . . ., and —j<m< j. The CS vector is obtained 
by displacing the vacuum with a SU{2)/U{1) coset element: 




with z — (i//|i/|)tan|i/| and normalization coefficient M — {1 -\- zz) K The CS manifold 
M = SU{2)/U{1) fa S^ CP^ is isomorphic to the 2-sphere S^, with 0<^<7r,0<<^<27r 
polar and azimuthal angles respectively. It is also parametrized as the complex projection 
plane CP^ with coordinate M. 3 z = tan(i^)e~*'^. 
The symbols of the generators are 

(J+) = 2j-^ , (J-) = 2j 



1 + zz 1 + zz 

Based on the general scheme the above generator symbols give rise to the line element 
ds'^ — (1 + zz)~'^dzdz for the complex tangent plane A4. The (over) completeness of the non- 
normalized CS invokes the Haar invariant measure on A4, dfi{z) = ^'^(1 + zz)^'^d'^z. The 
CS vector space is now the constant curvature {R = 1) sphere equipped with the canonical 
1-form 9 = ijj^dz + cc, which furnishes a canonical symplectic 2-form u = U (jf^fp-, and 

the associated Poisson bracket becomes {f,g} = ^(^+^^) [dzfdzQ — dzfdzg). 
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3 Generalized discrete self-trapping equation 



The classical Hamiltonian of the generalized discrete self-trapping (GDST) model reads 



/ / 
Hci. = E(^o|v4,f - ^|A,f ™) - £ XjkA.Ak , (11) 

where j = 1, 2, . . . , / counts the number of complex mode amplitudes Aj{t) and its complex 
conjugate Aj{t). The equation of motion 

/ 

is derived from the canonical Poisson brackets 

{Aj,Ak} = i6jk, {Aj, Ak} = {Aj,Ak} = (13) 
using the equation of motion Aj = {H, Aj}. 

The quantization of this canonical Lie-Poisson algebra proceeds with the correspondance 
rule ^ A — ^ b,A* — > h\ where 6^ and h are the canonical creation and annihilation 
operators respectively. To account for problems of ordering, well known in the naive quanti- 
zation rule, we shall here confine ourselves to two ordering rules, namely the normal ordering 
(NO) and the symmetric ordering (SO). SO is a classically motivated symmetrization of non- 
commuting operators and should be expected to give the right classical limit when — >■ 0. 
However here we study the quasiclassical limit of the quantum equations of motion by means 
of the CS symbol of operators and in this case the NO is the appropriate ordering. (For 
further discussion on the difference between semi- and quasiclassical ordering see refs. 89 
and 142 in |2|], ff p.64-65). 



For the SO quantization rule the correspondence |TB, ^ 

777 1 

l^r™ ^LU-2x) , (14) 

is valid for quantization of arbitrary positive powers of the modulus of the complex amplitude; 
here = : (b'^b)'^ : = b'^^b'' and Lm{-) is the Laguerre polynomial of zero order. Then, using 
the symmetric ordering quantization rule and after a constant shifting by ^uq — 7 ''™^"^''' , the 
Hamiltonian of the GDST model takes the form 

/ m—l 
j=l n=l 

Where /zM = (19^ f ^ 1 f , and 

/ / 
H^o = i2i^ob]bj - ^b^bf) - J2 \kh]bk , (16) 

is the Hamiltonian obtained by the normal ordering quantization rule. We observe that the 
symmetric ordering rule adds m — l additional terms in the Hamiltonian with respect to its 
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NO quantization. However since for m = 2 it is reduced to the number operator which can 
be absorbed in the H^q by simply redefining the coefficient of the corresponding number 
operator term, we see that in effect the symmetric ordering modiffes the Hamiltonian only 
for m > 2, i.e. only for higher order nonlinearity. 

Next we obtain the Heisenberg equations of motion for the boson operators by using the 
Hamiltonian of eq. ([TsD 

/ m-l _ 

tb, = ujob, - ^h^r'hj - y: \khu - 7( E ^^:r^nh^^ ■ (17) 

fc=l n=l 

Making use of the factorized state vector \P >= >, where bj\Pj >= Pj\Pj >,j = 

1, /, we evaluate the CS mean value of the last equation to obtain 

/ m-l 

4=^o/?i-7|/?,f'"-'/3,- E A,-fc/3fc-7E/^^^^l/^^f""'/^^-- (18) 

{j^k)k=l n=l 

By means of the canonical Poisson brackets this set of equations can be derived from the 
Hamiltonian 

n = (mip) = E(^oiAr - -\PkD -ji^kim - 7E EVir^i/^^r" (19) 

fc=l ^ k^l k=l n=l 

The above Hamiltonian and the norm = J2j=i IPjl"^ ^ire the constants of motion for the 
dynamics issued by eq. (JTSI). We observe that the first two sums in eq. (|TUp combine to give 
the coherent state symbol of the i^NO of eq. ([T6|) , while the remaining part corresponds to the 
symbol of the terms appearing in eq. (|T^) due to the SO. This implies that NO quantization 
induces in the CS manifold Ai an essentially classical time evolution while the quantization 
by SO, due to the additional terms, induces a dynamical evolution on Ad which departs from 
classical dynamics. This of course makes a difference in the dynamics of the observables of 
the system as well as in the wavefunction distributions. At this point we should remark 
that the assumption of factorization of the state vector into CSVs implies that any normal 
ordered function of boson operators from each mode has CS mean value expressed by the 
same function, i.e. 

(a| : f{a], aj) : |a) = f{aj , aj) . (20) 

Also for each oscillator mode the Hermitian combinations x = -^(a^ + a) and p = ;^(a^ — a), 
i. e. the position and momentum operators, we obtain the minimum uncertainty property of 
the boson CSV, 

AxAp = ^ , (21) 

for the uncertainties (Ax)^ = ((x — (x))^) and (Ap)^ = {{p — {p))^)- We shall recall this 
property shortly for the Q-function. 

To illustrate these differences in dynamics due to ordering we shall employ the quantum 
mechanical Q-distribution function |2^ and the Poissonian distribution of the boson excita- 



tions of a CSV in the number state basis [|T^. For the Q-distribution of a p density matrix 
we define 

Qj = -mpj\Pr) = -mm\\ (22) 

71 71 
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with j = 1, . . . , f . Here we have assumed a pure state for each mode i.e. pj = \j3j ){l3j\, and 
the reference CS 1/3^) is labelled by the complex coordinate jSr = Xr + iyr- The Q-function 
then takes the form of a displaced Gaussian function viz. 



Qj{Xr, Vr) = -e~^-r-ReP^?-(yr-lmP,f _ (23) 



The center of the Gaussian (Re/3j, Im/3j) is the expectation value of the position and mo- 
mentum operators in the CS basis, i.e. (Re/3j, Im/3j) = {{P\xj\P) , {P\pj\P)), for j = 1, . . . , /. 
However there is a nonzero quantum mechanical uncertainty around this point since Axj = 
Apj = albeit it is the minimum acceptable one (c./. eq. (pTD ). 

In Figs. l(a,b) the Gaussian bells of the Q-functions for the quintic (m = 3) GDST 
trimer (/ = 3) system are plotted at a certain instant of time for SO and NO respectively. 
Here and in the following we use periodic boundary conditions and choose Xjk = 1 for nearest 
neighbouring sites, and zero otherwise. The centers of the Gaussian functions of Fig. 1 are 
plotted in Fig. 2 for the SO (Fig. 2a,b) and the NO (Fig. 2 c,d) cases, respectively. We 
note that in the SO case the trajectory of the bell corresponding to the initially excited site 
will always be confined in a region far from the origin, while the bells corresponding to the 
other two sites will stay close to the origin. In our framework, this corresponds to the well 
known self-trapping of the classical system. It is obvious that the appearance of the extra 
terms in the SO-equation enhances the self-trapping compared to the NO-case. 
Next we consider the Poissonian distribution of each boson mode in the coherent state, 

P^ = m\n)\' = e-^'^^'^-^, (24) 



where n = 0, 1, . . . enumerates the number states and \(3jf = (^Pj\bjbj\(3j'^ is the expectation 
value of the number operator which gives the average number of quanta in the CS basis. 
An illustration of how the self-trapping reflects itself in the Poissonian distribution and how 
it gets enhanced in SO quantization, is shown in Fig. 3a and 3b for the SO and NO cases 
respectively. The system is here a rather large (/ = 21) quintic (m = 3) GDST-system, and 
initially one single site (j = 21) is excited with the total excitation number = 10. As 
can be seen, for the particular parameter values chosen here the main part of the excitation 
stays at the initially excited site in the SO-case, while it spreads more or less equally among 
the modes in the NO case. 

To end this section, we note that in the simple case of the GDST dimer (/ = 2) with cubic 
(m = 2) or quintic (m = 3) nonlinearity, it is possible to derive an analytic expression for the 
critical value,7cr, of the nonlinear coupling coefficient for the appearance of self-trapping in 
eq. ([I8|). This is done in the standard way by considering the variable r = — |/32p. For 
the ordinary GDST equation (|I2|) , it was found in ref. for the cubic case and in ref. |^ 
for the quintic case, that f could be expressed as a cubic polynomial in r, and solutions in 
terms of Jacobi elliptic functions were obtained. In the case of eq. (jl^), containing both 
cubic and quintic nonlinearity, we find the same equation for r as in ref. |2^, provided that 
we replace 7 with 'j^'^^ + N'y^^\ where 7*^^^ and 7^^-' are the coefficients of the cubic and quintic 
nonlinear terms, respectively. Thus, defining jcr in the usual way as the smallest value of 7 
for which \f5jof is always larger than A^/2, given that the excitation initially was localized 
on site jo, we find the self-trapping condition ^' 't^"^''^^ ^ ^ from the explicit solution for 
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r{t) just as in ref. ||22|. Using the explicit expressions for 7*^^) and 7*^^-' in (|T8D and \jk = 1, 



the following expressions for 7cr are obtained in the NO and SO cases, respectively: 

NO ^ A so ^ ^ (2^) 

These expressions are plotted in Fig. 4 as a function of the total number of lattice quanta A^. 
In the limit of large number of excitations (phonons) the quantum lattice model, according 
to Bohr's correspondance principle, will behave classically. In particular the critical value of 
7 for the onset of self-trapping, which as we have also seen from previous plots is different 
for SO and NO, is expected to become equal in the classical limit of the quantum lattice 
models. In this limit the ordering of the operators, a genuine quantum characteristic, should 



not be of importance. Indeed the analytic expressions (|25|) for the critical 7's of the GDST 
dimer show, that when the number of quanta becomes large, the phenomenon of DST occurs 
asymptotically for the same values i.e. 'j^^P ~ 7^*-*, independently of the ordering rule. 



4 Validity of the factorization ansatz and quantum cor- 
relations 

Let a composite quantum system be described by a general pure state \ip) living in H = 
®i=iHi, the tensor product of Hilbert spaces of its / subsystems. In general the form of 
such a generic vector I?/') can be of three different types, each one signifying a specific kind of 
correlation among the subsystems. For vectors of the product form, = ® ® . . . ® 
\4>)'^ , the subsystems are decorrelated. If \ip) is a convex combination of product states i.e. 
IV') = '^i \^)] ® 1^)1 ® • • • ® 10)/, with Aj = 1, the subsystems are said to be classically 



correlated ||29|. For inseparable or entangled systems, the third case, the state vector lip) is 
a general superposition of states from its subsystems with no additional properties, that are 
further characterized by e.g. violation of Bell inequalities this is the case when genuine 
quantum correlations are developed among the subsystems. 

Although the existence and the oddities of quantum nonlocal correlations and of entan- 
glement have been pointed out since the early days of Quantum Mechanics [|], and have 
remained an active subject of research since then it is only recently (due mainly to 
some striking applications of entanglement such as quantum computation, cryptography, 
teleportation etc; see |12| and refs. therein ), that measures that quantify the amount 
of correlations of a given state of a composite quantum system have been studied. In the 
above framework of quantum correlations it is clear that our factorization ansatz is equiva- 
lent to the assumption that no quantum correlations will be developed in the course of time 
evolution. Still, as we have seen, due to the nonlinear form of the quantum Hamiltonian the 
equations of motions derived from that approximate ansatz differ from the entirely classical 
ones. This difference captures some of the quantum features of the system, so we would like 
to have a figure of merit of our approximation. To this end we introduce a single index, 
the correlation index, that measures the distance of the true state of the model from the 
factorized one that is assumed, and so gauges the validity of our approximation. 

Recall that the CS-vectors form a basis for the Hilbert space of a single site, so for the / 
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sites of the boson chain model a general state vector can be expanded as 

P({a,}, {ai}) ®ti dfiiia^}) (26) 
with dfi{{ai}) = lli-id'^ai. By means of the relations ( [^ ), 

fiP) = I f{a)exp[-i{\af + \/3\')+M£a, (27) 
/(^) = f f{a){(3\a)d'a, 



which imply that the overlap of CS behaves as a delta function for the integrable 

functions of the Bargmann-Hilbert space, we can compute that {ip) is normalized provided 
that / \P{{ai}, {ai})fd^{{ai}) = 1. This shows that if P{{ai}, {ai}) = S{ai — (y{) ■ ■ ■ S{af — 
ap, then the general vector \ip) is reduced to the factorized state \ipref) = <S>f=i • 
We now introduce the correlation index e (c./eq. (|) ), 



e = V'i\l/J), \l/jref))=2-2m^ljref)\ (2J 
f - ^ 2 2- 

J dfi{{ai})P{{ai}, {ai}) J] exp[-i(|a[|^ + \aif) + 



2-2 



i=l 



Initially e(t = 0) = 0, since we set \ip{t = 0)) = \ipref)- In any subsequent time interval 
a < t < b, for which e{t = 0) ~ 0, our ansatz is justified. To determine e{t), we need to 
know the kernel P({aj}, {ai}). Substituting the expansion (|26|) in the Schrodinger equation 
we find that the kernel P satisfies the Fokker-Planck equation, 

thdtP{{ai}, {«,}) = n{ai, -da, + iai)P({«,}, {ai}) , (29) 

where Ti. is the analytic image of the Hamiltonian operator under the substitutions — *• 
ai , a] ^ —Bai ■ There are a number of approximate and analytic solution techniques avail- 
able for such equations [^, most accessible e.g. is the technique of continued fractions. 
However the problem of quantum correlations in nonlinear lattice models is interesting even 
in the dimer case, and as far as we know it has not been addressed. In such a case the 
conservation of the total number of excitation quanta in the lattice sites will make the num- 
ber of terms in a continued fraction method of solving the ensuing Fokker-Planck equation 
finite, since the dimensionality of the available Hilbert space and of the p-density matrix is 
determined by the number of quanta available. This makes possible an explicit calculation 
of the correlation index and allows to set a number of interesting questions concerning the 
relation of correlation/decorrelation versus selftrapping/non-selftrapping in the case of the 
quantum dimer (see e.g. for a full quantum treatment of the quantum dimer); however 
more work should be anticipated along these lines Finally, we note that for non-boson 



lattices, as the ones of the last appendix, the introduction of the correlation index can be 
done again based on the distance function defined for the respective CSV, c./eq. (||). 

5 Conclusions 

This work was inspired by the fact that for a linear quantum mechanical oscillator there 
is a coherent state vector defined in its Hilbert space, predicting that the motion will be 
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centered around the classical path with the minimum uncertainty. Here, this simple idea has 
been extended to some quantum nonlinear chain models. For quantum chain models with 
an oscillator attached at each site, the obvious assumption which has been followed here is 
that each site is in a CS ; the same idea however can be applied for multispin-like models 
with some advantages as we have shown. The essential geometric nature of our assumption 
has also been stated in terms of the symplectic and Riemannian geometry of the CSV 
spaces. Other types of special states, such as e.g. even/odd CSV (symmetric/antisymmetric 
combination of boson CS) which are associated with different geometries of their CS spaces, 
could have been considered; the final choice depends on the dynamical symmetry of the 
given Hamiltonian and on the specific form of the nonlinear interaction terms. 

Also, the treatment of the DST in the quantum regime as has been presented here 
shows that this is a phenomenon which continues to exist after quantization of the classical 
DST equations. However, as has been demonstrated by means of the CS method, its exact 
parameter dependence is crucially affected by the quantization scheme. Contrary to the 
NO case where the situation is essentially classical, in the SO rule DST gets even more 
pronounced from its classical form. We note that a similar statement can be made also 



concerning the Hartree approximation treated in [^, |18|, |3^, since this ansatz results into 



similar approximate dynamics as obtained here. This can be seen by comparing our eq. (18 



with eqs. (4.3) and (4.4) of ref. |T^. However, the Hartree ansatz is basically different from 
the factorization ansatz used here. Namely, geometrically it can be said that the two ansatze 
confine the dynamics in different subspaces of the total Hilbert space of the GDST model, 
moreover the Hartree wavefunction is not fully factorized. 

Finally, we should note that although the factorization hypothesis for the state vector has 
been effective in treating the dynamical evolution in the model-cases studied here, its validity 
should be questioned. This is a crucial point, since factoring the state vector of a compound 
quantum system implies a loss of the quantum correlations which may be developed during 
the course of the time evolution. The failure of a factorized state to account for such 
entanglement of quantum subsystems can however be quantified by means of a correlation 
index, as we have shown. In future work we plan to present a study of the entanglement 
phenomenon versus DST in a quantum dimer model. 
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Appendix 1: Modified discrete nonlinear Schrodinger 
equation 

The ability of the CS-method to extract the quantum mechanical characteristics of the time 
evolution for a given lattice Hamiltonian model is very much depending on the special form of 
the Hamiltonian itself. Occasionally the derived equations of motion will be reducible to the 
classical ones indifferently of the quantization ordering rule (at least for the common rules); 
let us take up the following case. The so called modified discrete nonlinear Schrodinger 
equation (MDNLS) |^ has been recently introduced as a model describing in the adiabatic 
limit polaron dynamics when the excitation is coupled to an acoustic chain of oscillators. 
Compared to ordinary DNLS dynamics of an excitation coupled to optical oscillators, the 
MDNLS contains some new formal features such as off-diagonal nonlinearity, and for a ring 
with /-sites reads 

iAj = V{Aj+i + Aj_i) - X(|A,+i|2 + + 2|A,f )A, . (30) 

Using canonical Poisson brackets, these equations can be derived from the Hamiltonian 

HcL = vj:{A,A,^^ + A,A,^,) - xj2i\A,^,\'\A,+2\' + |A,f )• (31) 

The quantization of the classical Hamiltonian for the MDNLS equation leads (after a constant 
shifting of the Hamiltonian) to the following expression 

/ 

Hso = H^o - J2ib]+2bH2 + b]+,b^+i + b]b,) , (32) 

where 

/ / 

b] b]) . (33) 

i=i i=i 
The first term in eq. ( |52| ) refers to the NO and the second one is obtained as before from the 
SO rule. Since the total Hamiltonian is written in terms of WH algebra elements defined at 
each site and by their products, we can follow exactly the same procedure as in the quantum 
DNLS case. Namely we assume the factorization of the state vector in terms of canonical 
CSVs. By the method of the previous section, and for the Hamiltonian Hso, this leads to 
derivation of a set of quasiclassical equations of motion satisfied by the CSV parameters Aj, 

%Aj = V{Aj+^ + Aj_i) - X(|A,+i|2 + + 2\A,\^)A, + ^XA, . (34) 

The above expression shows that the quasiclassical equation differs by a constant shift of 
the coefficients of the onsite variables from the corresponding classical equation. With the 
transformation Aj — >■ e~2*^*y4j they become the classical equations. The formal similarity of 
classical and quasiclassical equations allows to conclude the following: the particular form of 
the MDNLS Hamiltonian after quantization by NO and SO generates, (on the geometrical 
space singled out by the factorization of the model's quantum state in terms of boson CSV), 
an (almost) classical dynamics. This however as the previous sections show, should not be 
expected to be true in general. 
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Appendix 2: Fermionic polaron model (FP) 



To demonstrate the ability of the proposed method to study quantum dynamics for models 
others than quantized boson chains, we shall take up in this appendix the fermionic polaron 
(FP) model. It has been recently proposed to describe the interaction of electrons with 
optical phonons in one spatial dimension [|T5|. This same model has been studied from the 
point of view of QISM [|1|, ll] and the NSM 0, |6[. Under the assumptions made in ^ the 



Hamiltonian of the interacting electron-phonon system with periodic boundary conditions 
reads 



N N N 



H = a]aj - 51 ^(a]aj+i + aj+iflj) + V" ^ "■i^-i+i , (35) 
j=i j=i j=i 

where the fermion annihilation operator of a polaron at lattice site j, aj and its corresponding 
Hermitian conjugate Oj (the creation operator) generate the Grassmann algebra 

[ttj, ak]+ = [a], 4]+ = , [ttj, 4]+ = Sjk . (36) 

The occupation number operator is defined by rij = a^jaj, while g is a. parameter proportional 
to the overlapping integral and V stands for the electron-phonon coupling constant (see [13 
for details). Using the equivalence of the FP model to the XXZ model provided in |31| by 
means of a Jordan- Wigner transformation, we shall proceed to write the spin analogue of 
the fermionic model. The Grassmann algebra in eq. ( p6D is realized in terms of a tensor 
product of generators of the (spin 1/2 Pauli) algebra by means of the following Jordan- 
Wigner transformation 



i-i 

cj-exp I i7r^((T^ + ilj) 

k=l 

i-1 

k=l 



and its inverse 



= ^/^7 = (5li + ^|). (37) 

j-i i-i 
aj = ajexp{-i7i ^ Uk) , = 4exp(i7r ^ Uk) (38) 

k=l k=l 

a| = n, - 11, , (39) 

where af = ± zaj, and aj, aj, cr| are the Pauli matrices and Ij is the unit matrix at site 
j- 

Using eq. (|37|) the transformed Hamiltonian of eq. ( pS]) reads, 



where 



H = ^ie + ^) + {e + V)j:a] + H^^^, (40) 
^ ^ j=i 



H^xz = vj: a]a]^, - 2^7 E(^J^J+i + ^l^jVi) , (41) 
i=i i=i 
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or 

N N 

Hxxz = VY: ^-^Ui - 9T.i^t^J+i + ^7^/+i) • (42) 
i=i i=i 

We now proceed with the Hxxz psnt of the Hamiltonian as the other terms are constants 
of motion. Following same tactics as for the boson models we acknowledge the fact that 
the Hxxz Hamiltonian is embedded in the <^jLiSUj{2) algebra and make an appropriate 
ansatz about the form of the wavefunction. Namely, we assume it is written in the form 
\z) = <^f=i\zj >, where \zj > is the su{2) coherent state correspoding to the j-site of 
the chain. To determine the dynamics of the complex amplitudes z^s we use the relation 
z {j — Jo) = (</-), which is obtained from eq. (p^O]) . Then the time derivative of this relation 
for j = 1/2, gives the evolution equation of the complex amplitude at each site by means of 
the CS symbol of the Heisenberg equations, viz. 



z. 



d/dt (aj) (1/2 - a|) - {aj) d/dt (l/2 - a| 
1/2 - ^1 

aJ^H]) (1/2 -a^)^Uj)i\a^^^H_ 



,2 ■ (43) 

1/2 - a|) 

Straightforward evaluation of the CS symbols of the involved operators gives the evolution 
equation, 

1 



IZj 



l + |^,-_l|2)(l + |^,.+l|2) 

- Vzj{l - 2\zj-i\^\zj+i\^) + {2g - 3g\zj\^)[zj^i{l + \zj+i\^) + Zj+i{l + \zj^i\^)] 
+ gz][z,.i{l + \z,+i\^) + z,^i{l + \z,.i\^)]] (44) 

Returning to the original model we compute the CS mean value of the fermion operators by 
use of eqs. (|1^) and (pT]). This yields 

{z{t)\a^^\z{t)) = n(-2) {z,{t)\al\z,{t)) {z,{t)Wl\z,{t)) = C,(t)e-^^^W , (45) 

k=l 

where 



^'W-TTrrlrFnl-^E. (46) 



with Zj{t) = \zj{t)\e'^^^^\ 



Also, {z{t)\ak\z{t)) = (^z{t)\al.\z{t)j^ k = 1,...,A^, while the average value of the fermion 
number operator is given by 

{z{t)\n,\z{t)) = {z{t)\a^^\z{t)) + i = . (47) 

^ / 2 l + \zj{t)\ 

A final remark concerns the applicability of the CS method to the XXZ model; it has 
recently come to our attention that this method also can treat [^] the XXZ model of higher 
spins. 



13 



References 

[1] J. S. Bell, Physics (N.Y) 1,195 (1964). 

[2] F. A. Berezin, Comm. Math. Phys. 40, 153 (1975). 

[3] H. A. Bethe, Z. Phys. 71, 205 (1931); C. N. Yang and C. P. Yang, Phys. Rev. 150, 321 
(1966), and references therein; M. Gaudin, Phys. Rev. Lett. 26, 1301 (1971). 

[4] R. K. Bullough, N. M. Bogoliubov and G. D. Pang, in Future Directions of Nonlinear 
Dynamics in Physical and Biological Systems, ed. by P. L. Christiansen et al. (Plenum 
Press, New York, 1993), p. 217. 

[5] P. L. Christiansen, M. H. Hays, M. F. J0rgensen and A. C. Scott, Phys. Lett. A 201, 
407 (1995), and references therein. 

[6] P. A. M. Dirac, The Principles of Quantum Mechanics (Oxford University Press, Lon- 
don, 1947). 

[7] J. C. Eilbeck, P. S. Lomdahl and A. C. Scott, Physica D 16, 318 (1985). 

[8] A. Einstein , B. Podolsky and N. Rosen, Phys. Rev. 47,777 (1935) 

[9] D. Ellinas, J. Phys. A: Math. Gen. 26, L543 (1993); 

in Proc. "International Symposium on Coherent States: Past, Present and Future", Eds. 
D. H. Feng et. al (World Scientific, Singapore, 1994), p. 139; 

in Proc. "3rd Colloquium on Quantum groups and Physics", Czech. Jour. Phys. 44, 
1009 (1995). 

[10] D. Ellinas, Proc. "IV Workshop on Physics and Computation, PhysComp96" p. 
108. Eds, T. Toffoli et. al. New England Complex Systems Institute 1996, also in 
|http:/ /pm. bu.edu/PhysComp96/ . 

[11] D. Ellinas, M. Johansson and P. L. Christiansen, work in progress. 

[12] D. Ellinas and E. G. Floratos, J. Phys. A: Math. Gen. 32 L63-L69 (1999). 

[13] D. Ellinas, P. Maniadis, " g-Symmetries in DNLS-AL chains and exact solutions of 
quantum dimers" , |quant-ph/9907014|, submitted for publication. 

[14] V. Z. Enol'skii, M. Salerno, A. C. Scott and J. C. Eilbeck, Physica D 59, 1 (1992). 

[15] V. K. Fedyanin and V. Yushankai, Teor. Mat. Fiz. 35, 240 (1978); 
V. G. Makhankov and V. K. Fedyanin, Phys. Rep. 104, 1 (1984). 

[16] R. Gilmore, Lie Groups, Lie Algebras and Some of Their Applications (J. Wiley, New 
York, 1974). 

[17] R. J. Glauber, Phys. Rev. 130, 2529 (1963). 

[18] M. H. Hays and A. C. Scott, Phys. Lett. A 188, 21 (1994). 



14 



[19] S. Helgason, Differential Geometry, Lie Groups and Symmetric Spaces (Academic Press, 
New York, 1978). 

[20] T. F. Jordan and E. C. G. Sudarshan, Rev. Mod. Phys. 33, 515 (1962); ibzd 34, 377 
(1962); 

N. L. Balazs and B. K. Jennings, Phys. Rep. 104 (#6) 347 (1984); 

M. Hillery, R. F. O'Connell, M. O. Scully and E. P. Wigner, Phys. Rep. 106 (1984). 

[21] G. Kalosakas, G. P. Tsironis, and E. N. Economou, J. Phys.: Condens. Mat. 6, 7847 
(1994). 

[22] V. M. Kenkre and D. K. Campbell, Phys. Rev. B 34, 4959 (1986); V. M. Kenkre and 
G. P. Tsironis, Phys. Rev. B 35, 1473 (1987); G. P. Tsironis and V. M. Kenkre, Phys. 
Lett. A 127, 209 (1988). 

[23] T. W. B. Kibble, Comm. Math. Phys. 65, 189 (1979). 

[24] J. R. Klauder and B. -S. Skagerstam, Coherent States (World Scientific, Singapore, 
1985). 

[25] V. E. Korepin, N. M. Bogoliubov, and A. G. Izcrgin, Quantum, Inverse Scattering 
Method and Correlation Functions (Cambridge University Press, 1993). 

[26] W. H. Louisell, Quantum Statistical Properties of Radiation, (Wiley, New York, 1973). 

[27] V. G. Makhankov, M. Agiiero Granados, and A. V. Makhankov, J. Phys. A: Math. Gen. 
29, 3005 (1996). 

[28] A. M. Perelomov, Generalized Coherent States and Their Applications (Springer- Verlag, 
Berlin, 1985). 

[29] A. Peres, Quantum Theory : Concepts and Methods (Kluwer, Dordrecht 1993); 

G. Mahler and V. A. Weberruss, Quantum Networks, (Springer- Verlag Berlin 1995). 

[30] J. P. Provost and G. Vallee, Comm. Math. Phys. 76, 2890 (1980); 
J. Anandan and Y. Aharonov, Phys. Rev. Lett. 65, 1697 (1990). 

[31] F. C. Pu and B. H. Zhao, Phys. Lett. A 118, 77 (1986). 

[32] H. Risken, The Fokker-Planck Equation, (Springer Berhn, 1989, 2nd ed.). 

[33] E. Schrodinger, Naturwissenschaften 14, 664 (1926). 

[34] A. C. Scott and J. C. Eilbeck, Chem. Phys. Lett. 132, 23 (1986); Phys. Lett. A 119, 
60 (1986). 

[35] A. C. Scott and P. L. Christiansen, Phys. Scripta 42, 257 (1990). 

[36] A. C. Scott, J. C. Eilbeck and H. Gilh0j, Physica D 78, 194 (1994). 

[37] H. Weyl, The Theory of Groups and Quantum Mechanics (Dover Publications, New 
York, 1950). 



15 



[38] E. Wright, J. C. Eilbeck, M. H. Hays, P. D. Miller and A. C. Scott, Physica D 69, 18 
(1993). 

[39] W. M. Zhang, D. H. Feng and R. Gilmore, Rev. Mod. Phys. 62, 867 (1990). 



16 



Figure captions 

Fig. 1. The Q-functions for the quintic (m = 3) GDST-trimer (/ = 3) plotted 
at the time-instant t = 76.4 for the cases of symmetric (a) resp. normal (b) ordering of 
the boson operators. Dashed line corresponds to the initially excited site j = 2, while 
solid line corresponds to the sites j = I and j = 3. (They are equal for all times due to 
the symmetric initial condition.) In both cases, the total excitation number = 10, the 
nonlinear parameter 7 = 0.055, and the linear coupling coefficient Xjk is unity. 

Fig. 2. The trajectories of the centers of the Gaussian bells in Fig. 1 for times t < 260. 
The parameter-values are the same as in Fig. 1. Fig. 2(a,b) (Fig. 2(c,d)) corresponds to 
sites 1 resp. 2 for the SO (NO) case. 

Fig. 3. The Poissonian probability distribution from eq. (^) plotted as a function 
of excitation number n and site index j for the quintic GDST (m = 3) with / = 21 at the 
time-instant t = 50. Fig. 3(a) shows the case of SO, while (b) corresponds to NO. In both 
cases, the site j = 21 is initially excited, and the total excitation number is = 10. The 
value of the nonlinear parameter is 7 = 0.05, which gives self-trapping in the SO-equation, 
but not in the NO-equation. 

Fig. 4. The analytic expressions of 7cr eq. ( p5|) for a GDST quantum dimer with quintic 
nonlinearity {i.e. / = 2,m = 3) plotted as a function of the total excitation number for 
normal (NO) and symmetric (SO) ordering. The value of the linear coupling coefficients is 
set equal to 1 or to 1/2 if we consider periodic boundary conditions. We notice that the two 
7's become asymptotically equal for the moderate value of iV = 10. 
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